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Abstract 


A  formal  perturbation  theory  is  developed  for  Hamilton 
operators  with  piirely  discrete  or  with  mixed  (discrete  and  con- 
tinuoTOs)  spectra.  Formulas  for  the  resol'vemt  operator  are  given 
which  exhibitB  the  singularities  and   the  discontinuities  of  the 
matrix  elements  of  the  resolvent  in  the  representation  in  which 
the  unperturbed  Hamiltonian  is  diagonal.  These  results  are  ap- 
plied to  obtain  a  solution  of  the  time-dependent  SchrVdinger 
problem.  Two  kinds  of  asymptotic  expansions  of  this  solution 
are  then  given.  One  expansion  is  valid  for  large  values  of  the 
time  I  the  other  is  valid  for  values  of  the  time  sufficiently 
large  and  for  a  sufficiently  weak  perturbation. 
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1,  Introduction 

The  aim  of  this  paper  is  to  investigate  the  fonaal  theory  of  processes 
like  scattering  of  particles,  or  the  capture  and  the  emission  of  a  particle  by  a 
system.  We  are  concerned  especially  with  the  case  when  the  energy  of  the  incident 
particle  is  close  to  the  energy  with  which  the  particle  can  be  bound,  so  that 
resonance  phenomena  occur.  The  particle  co\ild  be  a  photon  and  the  system  an  atom, 
or  we  can  deal  with  nuclear  particles  and  nuclei.  The  mathematical  situation  in 
all  these  cases  is  characterized  by  an  unperturbed  Hamiltonian  having  a  mixed 
spectrum  partly  continuous  and  partly  discrete,  where  some  of  the  discrete  eigen- 
values lie  inside  the  continuous  spectjrum.  The  problem  is  to  evaluate  the  effect 
of  a  perturbation  on  such  ein  unperturbed  system  and  to  solve  the  corresponding 
time-dependent  Schrodinger  problem,  giving  also  asymptotic  formulas  for  large 

values  of  the  time. 

Problems  of  the  kind  treated  here  have  been  the  object  of  investigations 

by  several  authors.  Most  treatments  are  only  approximate  to  a  certain  order  in  the 
strength  of  the  interaction  (usiially  to  second  order).  This  approach  is  presented 
for  example  by  DiracLJ  in  his  discussion  of  collision  problems.  However,  the  very 
formulation  of  the  problem  in  these  papers  is  such  that  it  is  not  at  all  clear  how 
one  should  proceed  if  one  wants  to  calculate  higher  approximations.  Our  solution 
of  the  problem  exhibits  the  perturbed  resolvent  operator  in  terms  of  a  closed  set 
of  equations.  From  these  equations  an  iteration  procedure  immediately  follows, 
which  allows  calculation  to  any  order  of  accviracy  in  the  perturbation  parameter. 
A  perturbation  theory  for  the  resolvent  which  allows  calculations  to  all  orders  has 
also  been  developed  by  Heitler  and  his  collaborators  (see  [6]  where  other  litera- 
ture is  quoted).  Their  work,  however,  does  not  take  fully  into  account  all  the 
significant  consequences  due  to  the  presence  of  a  continuous  spectrum.  In  addition 
there  appear  to  be  some  difficvilties  with  their  use  of  the  concept  of  the  diagonal 
part  of  certain  operators.  In  our  treatment  the  use  of  this  concept  is  completely 
avoided.  Also,  we  give  expressions  for  all  the  matrix  elements  of  the  resolvent 
including  some  omitted  by  Heitler. 
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Perttirbations  of  continuous  spectra  have  been  studied  by  Friedrichs^ -*. 
One  of  the  res\ilt3  of  Firiedrichs'  paper  is  to  give  a  precise  mathematical  meaning 
to  the  kind  of  approximation  involved  in  the  formulas  given  by  Weisskopf  and  Wigner '-  -^ 
in  their  treatment  of  natural  line  breadth.  In  Section  6  we  shall  show  that  his  re- 
snlt  follows  quite  easily  from  the  properties  of  the  resolvent  operator. 

The  proofs  given  in  the  present  paper  are  not  mathematically  rigorous j 
rather  we  shall  follow  a  kind  of  heuristic  method.  Thus,  instead  of  stating  con- 
ditions to  be  satisfied  by  the  perturbation  and  then  proving  that  certain  conse- 
quences follow,  we  shall  simply  assume  that  the  pertvirbation  is  such  that  the 
necessary  operations  can  be  performed.  P'urther,  we  often  present  our  material  not  by 
treating  the  most  general  case,  but  jjistead  by  treating  in  detail  certain  special 
cases  where  the  important  features  are  clearly  exhibited. 

We  shall  solve  the  problems  mentioned  above  by  reducing  them  to  the 
problem  of  the  construction  of  the  resolvent  operator  for  the  perturbed  Hamilton- 
ian  H»  The  resolvent  operator  is  defined  as  the  operator 

(1)  G(X)  -  (X-H)-^  , 

where  X  is  a  complex  parameter.  It  is  well  known  that  the  spectral  theory  of  the 

operator  H  can  be  based  on  the  consideration  of  its  resolvent  •  Here  we  shall 

only  jiistify  in  a  heuristic  fashion  some  properties  of  the  resolventi  similarly, 

we  only  give  a  heuristic  justification  of  the  formula  (3)  below,  which  is  aised 

in  later  considerationso 

Let  us  call  P  the  projection  operator  corresponding  to  the  discrete 
m 

eigenvalue  X  of  H,  and  ?-.  the  projection  operator  corresponding  to  the  value  E 

The  inqDortance  of  the  resol-vent  operator  in  the  solution  of  problems  of  the 
kind  considered  here  has  been  emphasized  by  Sch*(5nberg  *- -•  *  L-'J ,  jn  Q,]  one 
can  also  find  several  mathematical  references* 
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of  the  continuoxis  spectrum.  We  can  then  write 

'^    dE 


(2)         a(x)  -  £]  3y-  +  '  ^^ 


the  integral  being  extended  to  the  continuovis  spectrxun  of  H.  It  is  evident  from 
(2)  that  the  matrix  elements  of  G(X)  between  two  normalized  states  are  analytic 
functions  of  X  which  exhibit  polar  singularities  at  the  eigenvalues  X  and  a 
discontinuity  across  the  continuous  spectnan. 

Analytic  functions  of  the  operator  H  can  be  constioujted  in  terras  of  the 
resolvent  by  a  formula  which  is  the  analog  of  Cauchy's  formula  for  analytic  func- 
tions of  a  conqjlex  variable.  Indeed  one  can  easily  convince  oneself,  using  (2),  that 


(3)  f(H)  '     ^     i       f(X)G(X) 


dX 


Here  the  function  f(X)  must  be  regular  in  all  points  of  the  spectrum  of  H  and  the 
integral  is  extended  to  a  closed  cxirve  C  (consisting  of  one  or  more  closed  loops) 
svich  that  all  eigenvalues  of  H  lie  to  the  left  of  C  and  all  singularities  of  f(X) 
lie  to  the  right.  Of  course,  if  the  spectrum  of  H  extends  to  infinity  in  the 
X-plane,  then  the  curve  C  must  also  extend  to  infinity,, 

We  shall  apply  formula  (3)  to  the  case  f(X)  ■  e~^   which  has  no  singu- 
larities in  the  whole  X-plane,  so  that  one  can  choose  any  contour  surrounding  the 
spectrum  of  H.  Formula  (3) then  gives  the  solution  of  the  tine-dependent  Schrodinger 
problem. 

Ptom  the  above  discussion  we  see  that  our  first  task  is  to  construct  the 
resolvent  of  the  perturbed  Hamiltonian  H  in  terms  of  the  unperttirbed  Hamiltonian 
and  the  perturbation.  This  is  done  in  Sections  2  and  3  for  a  p\irely  discrete  spec- 
trum of  simple  eigenvalues,  and  is  generalized  in  Section  k  for  the  case  of  a  mixed 


-i*  - 


(discrete  and  continuous)  spectrum,  which  is  more  important  for  our  purposes. 

We  then  use  oxir  results  to  obtain  the  asymptotic  formu.l.as  which  give  the 
solutions  of  the  time-dependent  Schrodinger  problem  for  large  values  of  the  time 
(Sections  5  and  6)» 

2,  Perturbation  of  discrete  spectra 

In  this  section  we  shall  give  the  formal  theory  of  perturbation  for  the 
resolvent  operator  in  the  case  of  purely  discrete  spectra.  fVom  this,  formulas 
for  the  eigenvalues  and  the  eigenstates  will  be  derived. 

We  consider  an  unperturbed  Harailtonian  H  which  we  assume  to  have  a 
ptirely  discrete  spectrum.  The  eigenvalues,  denoted  by  E  ,  will  be  assumed  to  be 
simple.  The  total  Hamiltonian  H  =  H  +  V  will  also  have  a  disci^te  spectirura  with 
simple  eigenvalues,  provided  the  pertiirbation  V  satisfies  some  conditions  which 
we  shall  not  investigate  here.  We  are  interested  in  the  matrix  elements 
/e  |g(X)|E  \      «  G  (X)  of  the  resolvant  G(X)  "  (X  -  H)   in  the  representation 
in  which  the  unperturbed  Hamiltonian  is  diagonal.  These  matrix  elements,  considered 
as  functions  of  the  complex  parameter  X,  have  polar  singularities  in  the  eigenvalues 
of  H»  We  consider  that  eigenvalue  X  of  H  which  tends  to  a  particular  eigenvalue  E 
of  E  as  the  perturbation  tends  to  zero.  We  would  likB  our  formulas  for  G__(^)  *o 
exhibit  clearly  the  polar  singularity  for  X  =  X.  . 

To  this  purpose  we  introduce  the  operator  |  '(X),  which  is  defined  in 
the  Ir-  representation  by  the  equation 

Clearly  ("^  (X)  depends  \:^on  the  value  n  of  the  index  which  is  omitted  in  the  sum 
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occurring  in  the  right-hand  side,  but  we  shall  not  indicate  this  dependence  ex- 
plicitly, since  we  want  to  concentrate  on  a  particular  value  of  n«   Since 
the  polar  singularity  for  X  "  E  does  not  appear  in  the  equation,  | 
can  be  ass'umed  to  be  regular  in  a  neighborhood  of  E  .  Actually  we  shall  see  in 


Section  3  that  |  '  is  regular  for  X  =  \  . 


We  can  now  give  the  formulas  for  the  matrix  elements  G  (X)  in  terms 

rs 

of  the  matrix  elements  |  '(X).  The  proof  of  these  formulas  will  be  given  later | 
in  this  section  we  shall  state  them  and  shall  derive  some  consequences  from  them* 
For  the  diagonal  element  r  »  s  -  n  one  has 


i$) 


x-E  -  n(x) 


n  'nn 
"Hie  other  matrix  elements  are  given  by 


(6) 


m 


1 


"^      x-E .  r(x) 


n  'nn 


(r  /  n)  , 


(7) 


G..(X) 


ns 


^•K-  ir„(« 


c 


ns 


X-E. 


n  'nn^^  s 

and  finally,  for  r  and  s  both  different  from  n, 

6 


(8)   <>rs(^)  '    T^*    T^ 


r 


^     C(x)  C(x) 
Q(x)  •     "^       ^^ 


'rs 


x-V  C(x) 

n  'nn 


(s  /  n)  , 


1 


These  formulas  show  that  the  polar  singularity  of  G«g(X)  in  the  vicinity  of  X  ■  B 


When  we  say  that  an  operator  is  regular,  we  mean  that  its  matrix  elements  in  the 


H  -representation  are  regular. 
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is  given  by  that  root  of 

(9)  X  -  E  +  r  (X) 

which  tends  to  E  when  the  perturbation  V  is  made  to  tend  to  zero*  This  root 
(which  is  a  simple  root,  from  our  assiunptions)  gives  the  eigen\-alue  \  we  were 
interested  in.  We  expect  it  to  be  real;  and  this  is  indeed  the  case,  as  can  be 

shown  using  the  reality  properties  of  |  (\)« 

Incidentally,  one  seas  very  easily  by  giving  r  and  s  the  appropriate  values 
that  (5),  (6)  and  (7)  are  contained  in  (8)»  Therefore  (8)  is  valid  in  general, 
without  restriction  on  r  and  s*  Hovrever,  if  r  or  s  or  both  equal  n,  the  simplified 
forms  (5),  (6)  and  (7)  are  more  useful,  because  they  do  not  show  the  denominator 
\  -  E  which  actually  does  not  correspond  to  a  singularity. 

To  fLnd  the  eigenvector  IV^  of   H  corresponding  to  X  we  need  the 
corresponding  projection  operator 

(10)  ^  "  ^ii      I   °^^^  ^^    * 

K 
where  the  integral  is  over  a  closed  curve  encircling  X_  t)ut  excluding  all  other 

eigenvalues  of  H.  (Formula  (10)  follows  immediately  from  (2)«)  One  obtains 

from  o\ir  formulas  for  0  (X)  the  expressions 


(11)  p   -   -i-   A    ±_-- dX  , 
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(ll4) 


n   r  n   s 


Obviously,  application  of  P  to  the  eigenvector  \^    ^    of  H  gives  \\y  $   apart 
from  an  arbitrary  normalization  factor.  To  obtain  comparatively  simple  formulas 
it  is  convenient  to  choose  |E  N  to  be  normalized,  <^E  |eN  -  1,  but  \\,^ 
such  that 

This  means  that  we  take 

pIe  N 
(16)  |\,>  =   2^-    . 

With  this  choice  of  normalization,  it  follows  from  (12)   that 

(17)       <K,^>  -  ^  [;;  (V  (-''"'• 

This,  together  with  (15),  gives  the  coefficient  of  the  development  of  |X  N  as  a 
linear  combination  of  the  normalized  eigenvectors  of  the  unperturbed  Hamiltonian. 

The  formulas  given  in  this  section  can  be  used  to  set  up  an  approximation 
scheme,  when  the  perturbation  V  is  sufficiently  small.  The  question  of  the  exact 
conditions  Tonder  wMch  this  is  possible  will  be  investigated  in  a  future  report. 
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3»  Perturbation  of  discrete  spectra  (contlnttation) 

We  shall  now  prove  the  fonirulas  (5)  to  (8)  for  the  resolvent.  For  con- 
cise notation  it  is  suitable  to  write  the  equations  of  this  section  without  re- 
ference to  a  particular  representation. 

We  introduce  the  projection  operator 

(18)  P".  |E„><EJ 

which  projects  on  the  eigenvector  of  H  corresponding  to  the  eigenvalue  E  .  The 
equation  which  defines  |  ' (X)  can  be  written  now  as 

(19)  PCX)  -  V  +  V(l  -  P°)  G°(X)  PCX)  , 

where  G(X)"(X-n)"  is  the  vinperturbed  resolvent. 

It  is  \iseful  to  find  different  forms  for  the  relation  satisfied  by  |  '• 
We  first  note  that 

(1  -  p°)G*'n  -  G(x  -  H°-  v)(i  -  p°)G°  r 
(20)  »  g[(i  -  p°)r  -  v(i  -  p^)G°r] 

-  g[(i -p°)r  -  r+ v]  -  G[y-p°r]   • 

Here  we  have  used  the  obvious  fact  that  P  commutes  w^th  G  ,  and  the  equation  (19) 
satisfied  by  |  '.  Substituting  (20)  into  (19),  we  obtain 

(21)  p  -  V  +  VG(y  -  pT"]  . 

Now  we  remember  that,  by  the  very  definition  of  G , 

(22)  GV   -  G(X  -  H°)  -  1 
and 

(23)  VG   -   (X-H**)G-1. 
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We  first  substitute  (22)  into  (21)  and  obtain 

(2W  n     -    VG(X-  Ii°-P°r)     . 

Then,  making  use  of  (23),  we  can  write 

(2$)      p  -  (X-  H°)G(X-  H°-  P°r)  -  (\-  H°-  p°n)     . 

If  we  write  (2$)  in  the  H  -  representation,  it  is  seen  to  contain  all 
the  X'elations  (5)  to  (8)  which  express  the  resolTent  in  terms  of  |  '•  To  see 
this  more  explicitly  without  use  of  a  particular  representation,  we  can  multiply 

(25)  to  the  left  by  P°G°,  which  gives 

(26)  P°  -  P°G(X-H°-P°r)  . 

This  formula  is  clearly  equivalent  to  (^  and  (?)• 

Now  we  observe  that  although  X  is  a  complex  parameter,  the  operators 
G,  G°  and  |  '  are  functions  of  X  'with  Herraitian  coefficients'.  This  means  that, 
for  instance, 

(27)  "T  (X)  -  ru)  . 

Therefore  if  (26)  holds  then  also 

(28)  P°   -  (X  -  H°-  nP°)  GP° 

holds.  To  prove  this,  one  need  only  take  the  Hermitian  conjugate  relation  of  (26), 
noting  that  X  is  an  arbitrary  complex  number.  Just  as  X  is.  Of  course  the  same 
observation  applies  to  any  of  the  other  relations  between  operators  we  are  dealing 
with« 

Using  (28),  we  can  write  (25)  as 

(29)  r  -  (X  -  H°)G(X  -  H°)  -  (X  -  H°)  -  PpVV   , 
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which,  when  solved  for  G,  gives 

(30)  G  -  G°+  G°  [r  +  rp°GP°r]  g°  . 

Using  (26),  we  find  that  this  last  formula  beccanes  identical  with  the  relation 
(8),  which  we  wanted  to  prove, 

Finally  we  want  to  give  a  relation  for  | '  i^ich  involves  the  resolvent 
G  in  a  more  symmetrical  form  than  (21)  and  (2k)   do.  We  obtain  it  by  substituting 
into  (21)  the  relation 

(31)  r  -  (X-  H°-  PP"*)  GV 

which  follows  from  (2ii)o     In  this  way  we  get 

(32)  n  ■    V  +  v|g  -  GP°(X  -  H°-   r')P°G]  V       . 

We  can  now  give  an  argument  to  the  effect  that  the  operator  f""  (X)  is 
regular  for  X  »  \  .  Writing  (32)  in  the  IT-  representation,  and  taking  into  ac- 
count (5),  we  obtain 


(33)        i;;-  ^«^^,^n  [«-¥"] 


"x. ' 


This  formula  gives  f"*  if  G  is  known;  it  can  be  considered  the  inverse  of  (8)» 

Now,  as  we  are  going  to  show  immediately,  the  combination  of  matrix  elements  of 

the  resolvent  which  appears  in  the  square  brackets  in  (33)  is  regular  for  X  »  \* 

Therefore,  under  certain  assumptions  on  the  matrix  elements  V  ,  the  elements 
'  rs 

I  '(X)  are  also  regular  for  X  ■  X.  • 

To  investigate  the  behavior  of  G.  j,  -  -g — Si   ait  X  -  X^,  we  first  re- 

nn 


call  that  _ 

(31.)  0,^  -  ^  '^ 


-  u  - 


we  have 

^    ^ht°hn      ^kn°ld 

(35) 

G 

15^    ^-\       ^-\ 

^nn 

Since  for  a  s\ifficiently  weak  perturbation  we  certaiiily  bSLve   c  c   /  0,  we  can 
say  that  the  right-hand  side  of  (35)  has  a  simple  pole  for  X  =  X  ,  and  that  the 
coefficient  is 

c  .  c  c  c  . 
fi<.\  nt  nn  nn  nJ?      — 

c   c 

nn  nn 

But  tliis  is  the  coefficient  of   «■■_■.■-   in  G.  ^  |  we  see  now  that  the  polar  singu- 

n 
larlties  of  the  terms  occurring  in  (33)  in  square  brackets  cancel  when  one  takes 

the  difference. 


li*  Pertiirbation  of  mixed  spectra 

The  work  of  the  preceding  sections  can  be  generalized  without  much  dif- 
ficulty to  the  case  of  an  unperturbed  Hamiltonian  ;*iich  has  a  purely  continuous 
spectinira  or  a  mixed  spectrum  (disci*ete  and  continuous).  In  this  section  we  shall 
treat  this  latter  case,  which  is  more  comprehensive.  Specifically,  we  consider 
an  unperturbed  Hamiltonian  describing  the  system  of  a  jjarticle  and  a  scatterer 
(cf«  Dirac'--',  Chapter  8),  lAen  there  exist  states  in  which  the  particle  is  bound 
to  the  scatterer.  The  unperturbed  Hamiltonian  will  then  have  a  continuous  spectinim, 
corresponding  to  the  situation  in  which  the  particle  has  a  given  momentum  and  the 
scatterer  is  in  a  given  state,  and  in  addition  it  will  also  have  a  disci^te  spectrum. 
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corresponding  to  the  situation  in  which  the  particle  is  boiind  to  the  scatterer. 
These  botind  states  will  be  assumed  to  be  orthogonal  to  each  other  and  to  the 
eigenstates  of  the  continu\un,  although  this  raay  only  be  approxiniately  true  for 
actual  physical  systems  (see  Dirac  •-  -',   ^1),  The  energies  of  these  bound  states 
may  lie  outside  the  continuous  spectrum,  or  may  correspond  to  values  of  the  con- 
tinuous spectrum.  In  the  latter  case  we  say  that  the  discrete  eigenvalues  are 
embedded  in  the  continuous  spectrum. 

We  now  introduce  a  perturbation.  The  effect  will  be  a  change  of  the 
discrete  states  and  their  energies;  the  continuous  spectnim  will  be  assumed  to 
remain  the  same  and  only  the  corresponding  eigenf unctions  vrill  xindergo  a  change. 
Formulas  for  the  perturbed  resolvent  which  exhibit  the  polar  singularity  corres- 
ponding to  one  particular  discrete  eigenvalue  can  easily  be  obtained  by  generaliz- 
ing the  formulas  given  in  the  preceding  sections  for  the  case  of  a  purely  discrete 
spectrum.  All  formulas  which  are  written  without  reference  to  the  PT-  representa- 
tion can  be  taken  over  as  they  stand.  The  essential  new  feature  occurring  now  is 
the  discontinuity  that  G°,  and  therefore  |  '  and  G,  exhibit  as  functions  of  X 
across  the  continuous  part  of  the  spectrum.  Let  us  consider  for  instance  the 
resolvent  G(X).  For  a  specific  value  of  the  energy  E  of  the  continuous  spectrum 
G(X)  takes  two  different  -values  according  as  one  approaches  the  point  E  from  the 
positive  or  the  negative  imaginary  side. 

An  important  consequence  of  this  discontinuity  is  that  one  can  try  to 
continue  analitically  the  matrix  elements  of  the  operators  involved  across  the 
continuous  spectrum  into  what  we  shall  call  the  second  sheet  of  the  plane  of  the 
complex  variable  X.  We  shall  assume  that  this  analytic  continuation  is  actually 
possible  to  some  extent.  In  a  certain  region  of  the  complex  X-plane  surrounding 
the  continuous  spectrum  we  then  obtain  two  branches  of  the  analytic  function 
G(X)j  using  conventional  notation,  we  write  G  (X)  for 
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that  branch  which  goes  into  the  resolvent  G(X)  of  the  first  sheet  on  the  positive 
imaginary  side  of  the  continuous  spectrum  and  g"(X)  for  the  one  which  agrees  with 
G(\)  on  the  negative  imaginary  side  of  the  continuous  spectrum*  The  sanie  notation 
Hill  be  used  for  \~^   and  G°. 

Turning  now  to  the  discrete  eigenvalues,  we  see  that  a  very  special  situa- 
tion can  occur  if  the  eigenvalue  E  is  embedded  in  the  continuous  spectrum.  By 
application  of  formulas  analogous  to  those  developed  for  the  case  of  discrete 
spectra,  we  wo\ild  expect  the  discrete  eigenvaliie  E  to  undergo  a  shift  of  position 
when  the  pertiu-bation  is  switched  on.  If  we  consider  E  as  a  pole  of  G  (X),  the 
perturbed  position  of  the  pole  will  be  given  by  that  root  of  the  equation 

(37)  X  =  E^+  <n|r*(X)|n> 

which  tends  to  E  when  the  perturbation  tends  to  zero.  It  will  be  shown  below 
that  in  general  this  root  has  a  negative  imaginary  part.  This  means  that  the 
perturbed  Hamiltonian  does  not  have  a  discrete  eigenvalue  at  all,  although  the 
matrix  elements  of  the  resolvant  have  a  polar  singularity  below  the  real  axis  in 
the  second  sheet.  The  discrete  eigenvalue  has  disappeared  under  the  influence  of 
the  pertxirbation.  This  situation  is  rather  typical  for  discrete  eigenvalues  era- 
bedded  in  the  continuous  spectrum  (see  Eriedrichs  ^ -"). 

In  the  case  jtist  discussed  the  matrix  elements  of  the  perturbed  resolvemt 
are  given  by  formulas  quite  analogous  to  (5),  (6),  (7)  and  (8).  We  write  them  here 
explicitly,  since  we  are  going  to  use  them  later.  For  siinplicity  we  consider  the 
case  when  there  is  only  one  simple  discrete  eigenvalue  and  this  is  embedded  in 
the  continuum.  We  denote  by  |e)>  the  eigenvector  of  H  corresponding  to  a  value  E  of 
the  continuous  spectiiunj  without  indicating  explicitly  the  other  variables  which  refer 
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to  "bhe  direction  of  the  motion  of  the  particle  and  to  the  state  of  the  scatterer* 
Thus  when  integrations  or  sums  with  respect  to  these  ■variables  occur  (as  for  in- 
stance in  formula  (li2)  below)  these  operations  will  also  not  be  explicitly  indicated. 
Furthermore,  we  simply  write  |n^  for  the  discrete  eigenstate  of  unpert'jrbed 
energy  E  .   We  then  have 

(38)  <n|G|n>    - 


X-E^.  <n|rU)|n>       ' 


(39)  <E|G|n> 


(ho)  <n|G|E> 


^      <E|r(X)|n> 


X-E-  <n|r(X)|n> 


X.E^-<n|r(X)|n> 
<n|rU)|E>     y^        , 


n 


m    <E|g|e'>    -     -5;-jr-     *     ^ 


<EiriE'>  + 


<E|r|n><n|r|E> 
X-E^-<n|r|n> 


X-E 


We  recall  here  also  the  definition  of  f"',  which,  written  explicit/  in 

the  Yi-   representation,  now  reads 

^   |E>  dE  <^E| 


(U2) 


ru) 


V  +  V 


X-E 


rex)   . 


The  integral  is  extended  to  the  continuous  spectrum  of  H  ,  and,  from  our  a8sviD?)tion 
that  only  one  discrete  eigenstate  is  present,  it  is  seen  to  be  equal  to  G  (X)(l-P  ;, 
where  P°-  |n><^n|« 

Now  we  wish  to  investigate  the  reality  properties  of  the  operator  fT'(X)» 
Of  course  we  assume  that  if   and  V  are  Hermitian  operators. 

We  recall  that  by  definition 

(i|3)  r(X)  -  V  +  V(l  -  P°)G°(X)  n(X)  . 
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Frora  this  we  can  prove 

(111*)  r(x)  -  ru)  , 

i.e.,  I*  is  a  function  of  X  with  Hermitian  coefficients.  Taking  X  =  E  to  be  real 
and  outside  of  the  continuous  spectrum,  one  concludes  that  the  operator  [~^  (E)  is 
Hermitian. 

We  know  however,  that  in  the  case  where  Yr   has  a  continuous  spectarum  and 
E  lies  in  it,  P  (E)  has  two  different  values  P  *(E)  and  ^"(E)  obtained  by  ap- 
proaching the  real  axis  from  the  positive  or  from  the  negative  imaginary  side  re- 
spectively. In  this  case  (14i)  can  be  used  to  show  tte.t 

(U5)  r*(E)  -  r"(E)  . 

To  find  out  more,  we  use  the  equation 

(1*6)  r(x)  -  V*  r(x)(i-p°)G°(x)v 

which,  in  accordance  with  the  observation  following  equation  (28),  is  also  satisfied 
by  I  '•  The  expression 

[l  *   r(X)(l  -  P^G°(X)]  v[l  *   (1  -  P°)G^x')r(x')] 

can  be  siiiQjlified  by  use  of  either  (1^3)  or  (lt6).  The  equality  relation  of  the 
two  resulting  expressions 

[i  +  r(x)(i  -  p°)G°(x)]r(x')  -  r(x)[i  ♦  (1  -  p°)GV)r(x')] 

can  be  written 

(a?)    r(x)-r(x')  -  r(x)(i -p'^[g°(x) -G°(x')]r(x')  . 

If  in  this  formula  we  take  X  »  E  4-  iC-,  x  »  E  -  ±0-    and  make  use  of  (Uli), 
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WB  obtain 

r  (E  +  i(^)  -   r(E  -  iC-)     -     nCE  +  iC-)  -   P(E  +  ±(^) 

-     n(E  +  i^)(l  -  P°)[q°(E  +  iO^)   -  G°(E  -  i(?)ln(E  +  ±<r)     . 

Now  we  let  C    tend  to  zero*  Remembering  that 

(U8)  lia       Ig°(E  +  iCT)  -  Q°(E  -  iC^)1     -    -  2ni  5(E  -  H°) 

we  have  _^ 

r*(E)  -  r*(E)  — 

a?)  21 -     -  «    r*(E)(l  -  P°)5(E  -  H°)r*(E)        . 

In  words:  the  anti-Herraitian  part  of  |    '   (E)   is  a  negative  definite  operator  (with 
E  lying  in  the  continuoxis  spectrum).     It  follows  in  particular,  since 
5(E  -  H°)    -   |E><E|,  that 

(50)  Im  <n|r*(E)|n>      «    -  n|  <n|  P^CE)  |E>  |      <    0     . 

We  shall  assume  now  that  <^n  !["'*(£)  |E^  /  0  for  all  values  E  of 
the  continuovis  spectrum.  The  imaginary  part  of  <^n|j  '  (E)  |n^  is  then  nega- 
tive for  all  those  values.  This  means  that  the  equation  (37)  cannot  have  a  real 
root  lying  in  the  continuous  spectirum.  On  the  other  hand,  it  cannot  have  a  complex 
root  in  the  first  sheet  of  the  X-plane,  since  this  would  correspond  to  a  complex 
eip;envalue  for  the  Hermitian  operator  H.  It  follows  then  that  if  (37)  has  a  root 
which  approaches  E  when  the  perturbation  tends  to  zero  and  if  the  pertixrbation  is 
sufficiently  small,  this  root  must  possess  a  negative  imaginary  part  and  lie  in 
the  second  sheet  of  the  \-plane. 
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5«  Asymptotic  behavior  of  the  solution  of  the  time-dependent  Schrodinger  problem 
We  proceed  now  to  investigate  the  asymptotic  properties  of  the  solution 
of  the  time-dependent  Schrodinger  equation  for  a  pertvirbed  Hamiltonian  of  the  kind 
considered  in  Section  k»     This  problem  can  be  reduced  to  the  evaluation  of  the 
resolvent  operator  by  means  of  the  forinxila 

(51)  e-^"*  «  ^    I  G(X)  e"^'^*  dX 


k  j 


2HI 

which  we  already  discussed  in  Section  1*  We  remember  that  the  contour  of  integra- 
tion should  encircle  the  vAiole  spectrum  of  the  Hamiltonian  H.  If  the  spectrum 
eoctends  to  infinity  the  contour  will  therefore  also  extend  to  infinity. 

Specifically  we  are  interested  in  the  matrix  elements  of  the  operator 
e"    in  the  representation  in  which  H  is  diagonal.  These  matrix  elements  have 
the  physical  meaning  of  transition  amplitudes  between  eigenstates  of  the  unper- 
turbed Hamiltonian,  from  the  time  0  to  the  time  t.  Since  we  have  in  formulas  (38) 
to  (Ul)  the  corresponding  matrix  elements  of  the  resolvent  G,  our  problem  is  solved 

by  (5l)  in  the  form  of  a  complex  contoTir  integral.  This  representation  is  veiy 

— iHt 
well  suited  for  investigating  the  behavior  of  t  he  matrix  elements  of  e     for 

large  t. 

We  shall  give  two  kinds  of  asyu^jtotic  formulas.  The  first  kind  will  be 
given  in  this  section  and  is  obtained  by  letting  the  time  t  tend  to  infinity  while 
the  other  parameters  remain  fixed.  The  second  asymptotic  formula,  which  will  be 
disciissed  later,  gives  the  behavior  of  the  matjrix  elements  for  times  sufficiently 
large  and  for  a  siifficiently  weak  perturbation.  In  addition  the  matrix  elements 
which  refer  to  a  transition  to  a  state  of  the  continuum  are  evaluated  for  unper- 
turbed energies  close  to  the  energy  of  the  discrete  \mpert\jrbed  eigenvalue.  The 
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precise  mathematical  meaning  of  these  statements  will  be  given  later*. 

We  begin  by  writing  explicitly  the  fommlas  which  give  the  matrix  ele- 
ments we  want  to  evaluate.  The  simplest  case  is  that  of  the  matrix  element 

-iXt 
X-E^-<n|r(X)|n> 

Since  we  want  to  investigate  the  behavior  of  the  integral  for  large  t,  it  is  con- 
venient to  change  the  contour  of  integration  so  that  \  has  everywhere  a  negative 
imaginary  part  on  it.  This  is  possible  only  by  moving  the  upper  half  of  the  con- 


(52)        <n|e-^"*In>   -  ^  (T  ^"^'^"_ -r-  dX  . 


tour  through  the  continuous  spectrum  into  the  second  sheet  so  that  it  lies  below 
the  real  axis.  We  see  in  this  way  that  <\n|e~   |n/>  tends  to  zero  as  t-^oo. 

This  result  is  rather  obvious  from  a  physical  point  of  view. 

— iHt 
We  proceed  now  quite  similarly  for  the  other  matrix  elements  of  e    . 


We  have 

(53)   <E|e-i"*|n>  -  ^ 


''  -iXt 
e 


^ni 


^^    <E!r(X)|n>  ^— /,_^^,^  dX. 


J 


X-E^-<n|r(X)|n> 

Moving  the  contour  as  before,  we  cross  this  time  a  polar  singularity  for  X  ■  E. 
Since  the  integral  over  the  new  contour  vanishes  asymptotically  for  large  values 
of  t,  we  can  write 

(5U)    lim  e^*<E|e-"^*|n>  -<E|n*(E)|n> 


t->oo  E  -  E^-  <n|n*(E)|n^ 

Notice  that  beca\ase  of  the  way  we  moved  the  contour  of  integration  the  |  '  branch 
of  the  function  |  '(X)  is  the  one  that  occurs  in  formvila  (5ii).  We  have  obtained 
in  this  way  the  final  value  of  the  transition  amplitude  from  the  discrete  eigen- 
state  |n^  into  the  eigenstate  |E^  of  the  continuum.  Remembering  that  the 
denominator  has  a  complex  zero  X  near  to  the  value  E  of  the  energy,  ve   can  say 
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that  the  transition  angjlitude  just  mentioned  must  have  a  sharp  maximum  for  E  =  real 
part  of  X  • 

Finally  we  are  interested  in  the  matrix  element 

m     <E!e-i"*|E'>  =  <E|E'>e-^*+  ^  (^  <EIr(X)  |e'>  -i-r  dX  , 
•where  we  have  set 

(.6)  <.nx,ih'>   .  <Eir(x,iK'>  .  <^in^><-iriE->    ^ 

X-E_^-<n|r|n> 

In  the  same  way  as  before,  we  obtain  asymptotically 

(57)  <E|e-i«*iE'>  -     <E|e'>   .  e"^*    <E|R*(E)|e'>     -^ 

E-E 

^    -^   <E|R*(e')|e'>     e-^'*       . 
£  -  £ 

Remembering  now  that 

it(E-E')  , 

(58)  lim ,-   -   in5(E  -  E  ) 

t  -i^oo    E-E 

we  finally  have 

(59)  llm       e^*<E|e"^"*|E'>     =     <E|e'>   +   f-iT  -  in6(£-E ')"]<£  |R*(E)  |E*>     . 
t  -^00  LE-E  J 

The  procedure  we  have  used  to  pass  from  (57)  to  (59)  is  familiar  from  the 
theory  of  scattering.  Formula  {^6)   shows  that  the  matrix  elements  of  R  exhibit 
a  resonance  for  E  ■  real  part  of  X  .  The  physical  interpretation  of  the  results 
of  this  section  is  well  known.  The  real  part  of  X  gives  the  displaced  position 
of  the  center  of  the  line  vdiich  corresponds  to  the  unperturbed  energy  E  ,  while 
the  absolute  value  of  the  imaginary  part  of  X  gives  the  width  of  the  line. 
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6»  Second  asymptotic  expansion 

In  the  preceding  section  we  have  given  expressions  for  various  matrix 
elements  in  the  limiting  case  t  ■  oo»  Now  we  shall  give  different  asymptotic 
expressions,  which  show,  in  some  cases,  how  those  limiting  values  are  approached 
as  t  increases* 

We  shall  treat  here  explicitly  the  case  in  which  all  assun^tions  of 
the  previous  section  conceiming  the  spectrum  are  satisfied,  and  in  addition  the 
perturbation  V  is  such  that  the  matrix  element  <^n|v|n^  as  well  as  all  elements 
<Ce|v|E  ^  vanish.  This  case  has  been  treated  by  Friedrichs*-  -^  and  the  refined 
asymptotic  expansion  mentioned  above  was  first  stated  in  a  precise  mathematical 
form  in  his  paper.  We  shall  show  that  it  follows  very  easily  from  our  formulas. 

Owing  to  the  particular  properties  of  the  perttirbation  V,  the  equation 
(I42)  for  the  operator  P*  (X)  can  be  solved  immediately.  At  this  point  it  is 
useful  to  write  V  ■  ttf,   where  e  is  a  parameter  measuring  the  strength  of  the 
perturbation.  As  one  can  see  very  easily,  the  matrix  elements  of  \       are  then 
given  by 

(60)  <E|r|n>.-  6<E|'ir|n>       , 

(61)  <n|riE>    -  6<n|l/|E>      , 

(62)  <E|r|E*>    -  0  , 

(63)  <n|  r  In> 


2 
e 


r 


J 

where  we  have  made  use  of  tte  conditions 
(6I4)  <Elt/|E'>  -     0 

(6^  <nll/|n>   -    0     . 


<ni?r|E>  dE<E(7rln> 
X^TE 
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These  restrictions  are  not  imposed  on  the  perturbation  only  for  the  sake  of  sim- 
plicity. In  actixal  physical  problems,  a  correct  choice  of  the  vinpertiirbed  Ifamil- 
tonian  will  ensiire  that  (65)  is  satisfied,  i.e.,  that  the  first-order  shift  of 
the  line  vanishes.  Hovrever,  equation  (6I4)  will  in  general  be  replaced  by  the 

statement  that  the  matrix  element  <Ce|v|E  ^  is  of  second  order  of  sroallness 

2  [2] 

(proportional  to  e  )  (cf .  Dirac   ,  ^  5l)  •  I't  is  easily  seen  that  this  less 

restrictive  condition  would  sxiffice  to  obtain  the  results  which  we  are  interested 

in  deriving. 

Before  we  proceed  f-urther  we  must  briefly  consider  the  equation 

(66)  X    -    E^+<n|r''(X)|n> 

which  we  have  assvuned  to  have  a  complex  root  X  below  the  real  axis*  Since  the 
fvinction  <^n|  |  *  (X)  \n^     is  assumed  to  be  analytic  in  X  in  a  region  containing 
the  continuous  spectnim,  and  is  analy+dc  in  6,  we  expect  the  root  of  (66)  to  be 
an  analytic  function  of  6.  As  shown  in  the  Appendix,  this  root  can  be  constructed 
as  a  power  series  in  e,  in  terms  of  -C^nll"*  (X)|n/>  and  its  derivatives  with 
respect  to  X,  calculated  for  X  »  E  .  These  quantities  can  be  expressed  as  improper 
integrals,  extended  to  the  continuous  spectnun,  involving  the  matrix  elements  of 
the  pertTorbation.  An  approximation  procedure  for  X  can  be  worked  out  in  this 
way  Wxhich  gives  consistent  results  to  all  orders  in  e.  Here  we  need  only  the 

It  is  usually  stated  that  the  power  series  expansion  in  s  for  the  shift  of  an 
eigenvalue  embedded  in  the  continuum  breaks  down,  because  if  one  tries  to  cal- 
culate the  coefficients  of  terms  containing  powers  of  e  higher  than  the  second, 
the  rasults  turn  out  to  be  infinite.  However,  we  show  in  the  Appendix  that  if 
this  calculation  is  properly  carried  out  one  obtains  a  shift  which  is  complex 
and  analytic  in  s. 
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rather  trivial  terms  of  second  order,  which  are  given  for  our  special  perturbation  by 
(67)  \  •=    E^+  6^0  -  ie^p  +   [e^] 

^     <n|V|E>dE<E|?/|n> 


with 
(68) 
and 
(69) 


P.V. 


J 


E     -     E 
n 


p     -    n<n|V|E^><E^|2/|n>     >0 


As  indicated  by  the  symbol  P»V«,the  integral  which  gives  a  is  to  be  evaluated  as 
a  principal  value  integral. 

We  can  now  proceed  to  the  asymptotic  evaluations.  We  begin  again  with 
the  matrix  element  <^n|e~   jn^  ,  and  vdsh  to  find  its  behavior  as  t  tends  to 

infinity  while  at  the  same  time  e  tends  to  zero  in  such  a  way  that  the  product 

2 

6  t  remains  finite.  We  can  still  use  (52)  and  move  the  upper  half  of  the  con- 
tour below  the  real  axis.  HovTever,  since  when  e  -»•  0,  X —»•  E  which  is  real,  it 
is  convenient  novr  to  move  the  contour  across  the  polar  singularity  of  the  inte- 
grand for  X  =  X  .  One  obtains  in  this  way 

<n|e-iH*|n> 


(70) 


-i)^t    ^   1 


n 


>'! 


where  the  curve  C  lies  as  indicated  in  the  figure  and  the  residue  in  X  gives 
the  factor 

r 


(71) 


1 


dX 


J 


^   X-E^-<n|r'(X)|n> 
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complex   X-plane 


C 
C 


The  continuous  spectrum, extending  over  the  real  axis  from  zero  to 
infinity,  is  a  line  of  discontinviity  for  the  f\inction  <^n|G(X)|n^. 
Dotted  lines  correspond  to  curves  lying  in  the  second  sheet. 

Now  we  note  that  when  e  -^  0,  then  X  — »E  and  r  — >-<^n|n^  ■  1.  Tne  contribution 
of  the  integreil  over  C  vanishes  in  the  limit.  As  a  consequence,  with  the  above 
specified  limiting  process,  we  have 


(72) 

or 

(73) 


iE„t 


14    ^n^  y    I  -iHti  \     -it 6  a  -  te  B 
lim  e  "  \n|«    \^/    "  © 


This  formula  gives  the  exponential  decay  of  the  probability  amplitude  for  t  he 
state  |n^  to  remain  unchanged.  We  want  to  emphasize  that  (73)  follows  from  the 
particular  way  in  which  we  assumed  the  perturbation  to  be  weak.  If  one  keeps 
the  strength  of  the  perturbation  fixed,  the  decay  of  a  state  wl.ll  not  in  general 
be  exponential. 
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To  obtain  the  asymptotic  evaluation  of  the  matrix  element  <^ E | e"   |nj/ 

we  shall  perform  a  limiting  process  in  which,  when  t  ->oo,  E  -»-E  in  such  a  way 

that  t(E  -  E  )  remains  finite.  In  addition  e  -♦  0  as  before,  so  that  te  remains 

finite.  The  asymptotic  formula  which  we  obtain  in  this  way  (formula  (76)  below) 

will  therefore  give  the  matrix  element  for  sufficiently  large  values  of  t,  for  a 

sufficiently  weak  perturbation  and  for  values  of  the  energy  sufficiently  near  to 

E  •  We  use  (53)  and  move  the  contour  as  before.  Then 
n 


lEt 
<E|e-^"*|l,>  -  <E|r*(E)ln> 


E-E_j-  <n|P*(E)|n> 


*<^irX'l">   V^  'n*  53 


c 


J 

With  the  specified  limiting  process  we  have 

lim  t(E  -  \j)     -    t(E  -  E^-  e^a  +  is^)     =     lim  t(E-E^-  <nl  n*(E)  |n>  ) 
and 

limt^/2  ^E|r*|n>     -    t^^6<E^[V|n>         . 

Therefore  we  can  state  that 

(75)       Xi.t-A,i^t<.|e-iHt|.>     .     .1A,<KJ.|„>     ^-«P["'^-Vf°-^'^P'] 

^  ^   n        /        t  E-E^-  6  a  +  ie'^pl 

or  _ 

e3q)[-iEt]  -  exp  -it(E  +  6^a)  -  te^p 

(76)<E|e-i«*|n>^   e<EjUln>      — ^--1 \ ^-5- ^  • 

^  Nn/  E-E-6^a  +  ie*^ 

n  ^ 

Finally,  a  similar  method  can  be  applied  to  the  evaluation  of  s^EJe"   |E  ^» 
Here  we  must  also  let  E  ~>"E  in  such  a  way  that  t(E  -  E  )  remains  finite.  We  do 
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not  reproduce  the  complete  calculation^  which  follows  the  same  pattern  as  for  the 
matrix  elements  evaluated  above.  The  result  is 

li^l  [e^XE|e-^«*|E'>-<E|E'>]  -  te2<Eji;  |n><n|ir  |E^ 

exp[it(E-.E  -  e^a  +  i6^p)j 


(77)   X 


r — ^ 

t(E-E^-  6  o  +  it 


or 


.6'p)t(E-E  ) 

exp[it(E  -  e')J 
t(E  -  E  -  e^a  +  ie^p)t(E  -  e') 


<E|e-"t|E'>~<E|E'>  .-«*  .     s^  <EjV|n>  <n|V|E^> 


t(E  -E  -  t'^a   +  i6^p)t(E'-  E  -  e'^a+ie'^p ) 


(78)   X 


exp 


JT.iEt] 


expl-it(E^+  6  a)-ts^p 


] 


(E-E  -  e^a  +  i6^p)(E-E')      (E-E  -  e^a  +  ie^p)(E'-  E  -  s^a  +  ie^) 
n         "^  n        "^      n        "^ 

exp[-.iE  t] 


(e'-  E  -  6^a  +  ie^p)(E  -  E*) 


It  is  interesting  to  notice  that  if  in  the  asymptotic  formulas  (76)  and  (78)  we 
now  let  t  tend  to  infinity,  there  resxilt  expressions  similar  to  those  given  in 
(5I4)  and  (59)  respectively,  only  with  E  substituted  by  E  everywhere  except  in 
the  first  term  of  the  resonance  denominator.  Also,  if  in  (73)  we  let  t  tend  to 
infinity,  the  expression  vanishes,  just  as  the  exact  matrix  element  given  in  (52) 
does* 
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Appendix.  Solution  of  equation  (66):  X  »  E  +  i^^nl  [~^  (X)|n^ 

In  the  case  of  a  discrete  eigenvalue  embedded  in  the  continuous  spectrum, 
the  perturbed  resolvent  has  a  pole  for  X=X,  withXa  solution  of  (66).  In 
this  Appendix  we  shall  investigate,  somewhat  more  thoroughly  than  in  the  main  text, 
the  problem  of  solving  (66)  and  in  particular  the  construction  of  X  by  an  expan- 
sion in  powers  of  the  strength  of  the  interaction.  We  shall  concentrate  on  the 
simplified  case  treated  in  Section  6,  so  that  (66)  takes  the  fom 

f  \<n\is\Eyf 


X  »  E^+  6^ 
n 


J 


X^TI ^ 


We  are  led  in  this  way  to  the  investigation  of  the  solutions  of  an  equation  of  the  form 


(A«l)  X  -  Xq+  o 


r      |f(E)|2 


T^rr   ^ 


J 


where  X  lies  inside  the  inteinral  [a,b3  •  In  actual  physical  problems  the  continuous 
spectrum  extends  from  zero  to  infinity,  but  for  reasons  of  simplicity  we  shall  con- 
sider here  the  case  of  a  finite  interval.  We  also  assuiae  that  the  function  tiX^ 
is  analytic  and  regular  in  a  region  R  of  the  complex  ^-plane  containing  the  in- 
terval [a,b^  in  its  interior.  The  function  7(^)  is  then  an  analytic  function  of 
X,,   let  us  call  it  g(^).  Clearly  g(^)  is  regular  in  a  region  S,  \diich  can  be  ob- 
tained from  R  by  a  reflection  with  respect  to  the  real  axis.  The  product 

(A.2)  f(^)g(^)  -  h(d 

is  regular  in  the  intersection  T  of  R  and  S,  which  is  a  region  symmetric  with 
respect  to  the  real  axis.  For  ^  »  E  real,  we  have 

(A.3)  h(E)  .  f(E)g(E)  -  f(E)f(ET  -  if(E)|2  . 
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We  can  write  equation  (A,l)  in  the  form 

b 

(A.U)  ^  -    V«   f  1^  ^ 

and  we  know  now  that  the  function  h(^)  is  regular  in  a  region  T  of  the  con^Jlex 
^-plane  which  contains  in  its  interior  the  interval  jajb]  of  the  real  axis. 
Consider  now  the  function  of  X 

(A.5)  F(X)  -  j   |i£^  dE   . 

It  is  analytic  in  X  and  its  representation  (A. 5)  shows  that  it  is  regular  in  the 
whole  X-plane  with  the  possible  exception  of  the  interval  [a,b].  Actually  it  is 
easily  seen  that  this  interval  is  a  line  of  discontinuity,  the  Jump  across  the 
real  axis  at  the  point  E  being  given  by 

(A.6)     lim    TfCE  +  iCT)  -  F(E  -  i<^)]  -  2Hi  h(E)  . 

(r->  0 

As  a  consequence  of  the  analyticity  properties  of  h(^),  one  c?m  continue  r(x) 
across  the  discontinuity  into  a  second  sheet  for  the  complex  variable  X*  For 
instance,  if  one  approaches  the  discontinuity  interval  from  the  positive  imaginaiy 
side,  one  can  evaluate  the  integral  along  a  curve  C  connecting  a  and  b  and  lying 
below  the  real  axis  but  inside  the  region  T.  In  this  way  one  sees  that  F(X)  can 
be  continued  in  the  negative  imaginary  direction.  The  boiindary  of  T  poses  a  natural 
limitation  to  the  continuation  process*  The  endpoints  a  and  b  of  the  continuous 
spectrum  appear  as  branch  points  for  the  function  F(X)«  Following  a  notation 
already  employed  in  the  main  text,  we  shall  denote  by  F  (X)  the  branch  just  con- 
structed of  the  analytic  function  F(X),  which  is  regular  in  a  region  containing 
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the  open  interval  (a,b).  Similarly,  we  would  obtain  F"(X)  approaching  the  line 
of  discontinuity  frcan  the  negative  imaginary  side. 

Equation  (A.l;)  can  be  written  now  more  precisely  as 

(A.7)  X  =  X  +  a  F*(X), 


where  the  function  F  (X)  is  regular  in  a  region  containing  X  in  its  interior. 
It  is  clear  that  (A,7)  admits  a  solution  X  =•  X  which  is  analytic  in  a  and  ap- 
proaches X  as  a  approaches  zero.  This  solution  can  easily  be  given  as  a  power 
series  in  a.  Successive  differentiation  of  (A.7)  gives 


(A.8) 


l^X 


da' 


dF*  dX 


dX   da 
a 


iV 


+  a 


dx: 


dP*  d^X 


+  a 


a 


dX   da 
a 


and  so  on.  Setting  a  =  0,  we  get 


>  -  0 


«•  X 


(A.9) 


dX 

0 

da 


F*(X^) 


d^X. 


dfl 


a  •=  0 


2  ^* 


?*(X^) 


and  so  on.  Thus,  to  obtain  the  power  series  expansion  for  X  ,  one  only  needs  to 
know  the  values  of  F  (X)  and  its  derivatives  with  respect  to  X,  calculated  for 


X  =  X  •  Remembering  (A. 5)  and  (A.3),  we  can  write  these  qviantities  as 
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(A.1C) 


d°F^X) 


dX" 


b 

r 


x=x 


G-  -^  0+ 


|f(E)r  dE 


J 


^   (X  +  iO-  -  E) 
ft   o 


■:^        (n  -  0,1,2,...) 


It  is  important  to  notice  that  the  limits  occurring  in  the  right-hand  side  of 
(A.IO)  actually  exist,  although  X  lies  inside  the  interval  of  integration. 
Indeed,  since  h(^)  is  regular  in  the  region  T,  one  can  integrate  along  the  con- 
tour C  mentioned  above  and  obtain 


b 
r 


(A.ll) 


lim 
0-  ->  0+ 


|f(E)r  dE 


r 


a  <  V  ^^  -  ^^ 


i?T 


J 


hCc)  d^ 

0  V^>^ 


Concerning  reality  properties,  we  notice  that  all  quantities  given  in 
1.  This  one  can  see  from  the  fact  that  the  complex  cc 
woxild  be  obtained  by  a  formula  similar  to  (A. 10)  only 


(A.IO)  for  n  >  0  are  real.  This  one  can  see  from  the  fact  that  the  complex  con- 

dV 

jugate  of 


dX 


n 


X-X. 


with  0"  tending  to  zero  through  negative  values.  In  the  analog  of  (A.ll)  one 

t 

would  then  have  an  integral  over  a  ciirve  C  lying  above  the  real  axis.  However, 

this  curve  can  be  moved  to  coincide  with  C,  since  the  residue  at  the  pole  X 


vanishes  for  n  >  0.  Hence 


d"F* 


dX 


n 


coincides  with  its  complex  conjugate. 


X«X. 


For  n  =  0,  one  has  the  familiar  expression 


(A.12) 


F*(X^)  -  P.V. 


"^      !f(E)|^dE 


mr-  -  i"if(VI  > 


J 


a 


where  the  integral  has  to  be  evaluated  as  a  principal  value  integral,  "a  this 
case  the  negative  imaginary  part  -n|f(X  )|  results. 

It  is  clear  frcan  the  above  work  that  the  case  in  which  \    lies  outside 

o 

of  the  interval  Ia,b]  would  lead  to  the  same  formulas  (A. 9)  and  (Adfl)  for  the 
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coefficients  of  the  power  series  expansion  of  \  .  However,  in  this  case  one  can 
evalT;iate  the  limits  in  (A,  10)  by  j\ist  setting  C"  •»  0  in  the  integrand.  The  re- 
sulting f orm\ilas  are  then  identical  vdth  s  tandard  perturbation  theory*  This 
sho>7s  that  the  case  treated  here  of  a  discrete  eigenvalue  embedded  in  the  con- 
tinuous spectrum  can  be  considered  as  solved  essentially  by  standard  pert\irbation 
theory  if  one  adds  the  provision  that  the  singular  integrals  occxirring  must  be 
evaluated  by  the  limiting  process  specified  in  (A. 10) •  This  limiting  process 
extracts,  so  to  speak,  the  finite  part  of  integrals  of  the  kind 


f       |f(E)|2dE 


J 


„  (^  -  A) 
a   o 


with  X  contained  inside  the  interval  ja,b"j»  It  wovild  be  easy  to  see  that  the 
present  definition  of  the  finite  part  of  an  integral  is  identical  with  the  de- 
finition given  by  Hadaaard  I-  J  and  widely  used  in  the  mathBmatical  literature. 
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